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1. Introduction 

The Wilson line is a requisite to elucidate the high-energy collision of partons in the eikonal 
approximation H. Especially in case of scattering between a light projectile and a heavy 
target such as the electron-hadron collision, proton-nucleus collision, deuteron-nucleus col- 
lision, and also nucleus-nucleus collision in the forward or backward rapidity region due to 
small-x evolution, etc, the scattering amplitude is expressed in terms of the Wilson line 
representing partons which reside in a light projectile and travel through random color 
fields from a heavy target ||, ||| . It is not the individual particles inside the target but 
its surrounding fields that the projectile can probe. Such a description is analogous to the 
Weizsacker- Williams approximation in which electron is viewed as equivalent photon. The 
notion of non-Abelian analogue of the Weizsacker- Williams field has been well developed, 
which is called the Color Glass Condensate (CGC) ||, |j in the field of high-energy QCD. 

The McLerran-Venugopalan (MV) model assumes that the color charge density p is 
static (x + -independent), is a function of the transverse x±_ and longitudinal x~ coordi- 
nates, and distributes randomly at each spatial point. Its magnitude squared \p\ 2 should 
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be proportional to the transverse density of partons consisting of a heavy target, which 
is commonly denoted by /i 2 in the model. One can compute the scattering amplitude 
by taking the Gaussian average of Wilson lines embodying the projectile given a cer- 
tain n relevant to the experimental condition. Since the explicit form of the non-Abelian 
Weizsacker- Williams field is known ||, the above mentioned is a doable calculation. 

In fact, one can find evaluation of the Gaussian average of Wilson lines in literatures ]3|, 
H> 0) H) H> m different contexts and thus with different color structure, representation, etc. 
We here aim to derive more general formulae, which provides us with useful implements 
to describe high-energy collisions. The most general form is, as easily anticipated, too 
complicated to handle directly once the number of Wilson lines is more than four, as we 
will encounter later in this paper. In that case we will attempt to simplify the expression 
under the limit of large N c where N c is the number of colors. We will see that a picture of 
the color dipoles instead of gluons naturally arises in the large- iV c limit. 

Just for clarity of what we will address, we prefer to use the terminology, "scattering 
amplitude" to signify the Wilson line correlator. That quantity is, however, not limited 
only to the scattering process but would appear in the process of particle production from 



the CGC background || [lC], ||, 11 1. Also, we would mention that the Gaussian average 



is not only limited to the MV model but is widely relevant to the CGC formalism with a 
Gaussian approximation fl2|| . Therefore, we believe that the potential application of our 
formalism should be ubiquitous in high-energy QCD. 



2. Gaussian average of Wilson lines 

Our goal is to derive the general expression of the Gaussian average or correlation function 
in terms of Wilson lines under random distribution of color source. In a physical terminol- 
ogy the correlation function represents the scattering amplitude of a bunch of particles and 
antiparticles traveling through random color source in the eikonal approximation. That is, 
the specific quantity of our interest in this paper is 

(U(x 1± ) l3iai U(x 2± )p 2a2 ■ ■ -U(x nA _) 0nan ) , (2.1) 

where the Greek indices are with respect to color in a certain representation r of the SU(iV c ) 
group. In the MV model the Wilson line is written by the non-Abelian Weizsacker- Williams 
field given as a solution of the classical Yang-Mills equation of motion; 



U(x±) = "Pexp 



/+oo 
dx~d 2 z_i_ G (x±-z±) p a (x~, z±) t a 
-oo 



(2.2) 



where t a 's are color matrices of the SU(A r c ) algebra in the r representation. We denote the 
time ordering operator in the x~ direction by V and the two-dimensional propagator by 
Go(xj_) which satisfies the Poisson equation, 

1)1 G (x ± ) = 5^(x ± ). (2.3) 



dx 



2 
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[U(x 1± )f3 iai U {x 2i _)(3 2a2 ■ U (x nA _)p nan ) 



Figure 1: Graphical representation of the correlation function in terms of Wilson lines which 
corresponds to the scattering amplitude of n-particles traveling through a heavy target with the 
random and dense color distribution. Each Wilson line stands for a parton which starts with the 
color orientation oti before scattering and ends up with [3i with a fixed transverse position x,_l due 
to the eikonal approximation. 



Figure [l] is the schematic picture of the average Q2,l| ) with color indices. The blob part 
is the target which provides the random and dense p a (x~,x±) from the target, where x~ 
and x±_ indicate the spatial point on the transverse (impact-parameter) plane and the 
longitudinal extent of the target respectively. It should be noted that x~ is regarded as a 
time variable for the projectile. Thus, the Wilson line ( |2.2j ) encodes projectile's multiple 
scattering off the CGC along the temporal x~ direction. 

It is assumed in the MV model that the average (• • • ) is accompanied by the Gaus- 
sian weight in terms of p a (x~ , £Cj_), whose dispersion specifies the typical model scale p in 
the standard convention, or in other words, the saturation scale Q s related to p up to a 



logarithmic factor (see Eq. (2.11) for our definition without logarithm) universally charac- 
terizes the hadron wavefunction. As we mentioned before, we will develop our method for 
the MV model for example, but the technique is applicable to any CGC calculation with 
a Gaussian weight function as adopted in Ref. |l! 
The explicit form of the Gaussian weight is 



u)(p) = exp 



2p\x- 



(2.4) 



In fact, the random walk in SU(iV c ) group space leads to the quadratic term |L3|] in the 
weight function, and besides, the cubic term |Q] which is sensitive to Odderon exchange 
but is beyond our current scope. 

The only necessary ingredient for our calculation in what follows IS 5 ELS Si matter of 
fact, the two-point function of p a which is spatially uncorrelated as 

(p a (x~,x±)p b (y~,y ± )) = 5 ab 5(x~ - y~)5 ( - 2) (x ± - y ± ) p?(x~) , (2.5) 

which contains the equivalent information as the weight function (|2,4|). 

Now we have finished the setup of the MV model, that is, we have explained the 
notation and the model definition in a self-contained manner. In the subsequent discussions 
we will proceed toward the general expression in the n-particle case step-by-step starting 
with the simplest case of one particle. 

2.1 One-point function; (U(x±)p a ) 

We aim to make clear our notation (which is the same as Ref. ||) first in a warming-up 
exercise though the average of one-point function is not physically relevant. Our treatment 
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Figure 2: Contraction of four sources; (a) tadpole type, (b) nesting one, and (c) overlapping one. 
Two dots connected by the wavy line are contracted to the same point (time) and only tadpole-type 
diagrams remain finite. 

and convention are parallel to Ref. ||. Here we introduce the Wilson line integrated over 
a finite range denned by 



U(b ,a \x±) = V exp 



-15 



/ dz d 2 z ± G Q (x ± -z ± ) p a (z ,z ± )t a 

J a~ 



(2.6) 



The limit of a - — > —00 and b~ — > +00 renders the above the Wilson line as defined in 
Eq. ( |2.2| ). We will expand the finite ranged Wilson line and compute its Gaussian average 
using Eq. fl2.5|) . The Taylor expansion of time-ordered exponential function leads to 

00 n „{,- 



(U(b-,a-\x ± )) = ]T(-uj 2 r f f[d 2 z i± G (x ± -z iX ) [ dsj" f" darj • ■ ■ f 

n=0 ^ i=l • /a Ja Ja ~ 



X <Pax fa , Zl±)pa 2 fa > *2_0 ' • • Pa„ fa > ^ ^ 



dz„ x 



(2.7) 



Here we can decompose (• • • ) into all possible contractions in case of the Gaussian average. 
Then, only the adjacent contraction making the tadpole diagram as shown in Fig. || (a) 
survives and other contractions as in Figs. || (b) and (c) vanish because of the delta- 
function in Eq. (|2.5| ) . Since the tadpole contribution is to be factorized as {paiPa 2 ' ' ' Pa n ) — 
(pa 1 Pa 2 )(Pa 3 ■ ■ ■ Pa n ) , we can rewrite Eq. into a form of the integral equation; 

(U(b~,a~\x±)p a ) = 

= Spa + Hs 2 ) 2 / dz^d 2 z 1± I dz^d 2 Z2±Go(x ± -z 1± )G (x ± -z 2± ) x 

X (Pai fa~ , Z 1± ) Pa2 fa , Z 2± ) > (t* 1 1' 12 ) ^ (U(z 2 , a~ \ X ± ) la ) 

= Sp a - —C 2 {r)8p 1 J d 2 z_ L Gl(x ± -z ± ) J dz~ p 2 (z~)(U(z~ , a~\x±) ia ) , (2.8) 

where we note that we used J 1 dz^ <5(fa — z^") = \ and (t a t a )p a = C 2 (r)5p a with the 
second-order Casimir invariant C 2 {r) in the r representation. We show the diagrammatic 
representation of this integral equation as 



+ -■ 




with the tadpole attached at {z ,x±). We can easily find the solution, that is given by 



(U(b ,a \x±)p a ) = U(b ,a \x±) 6p a = exp 



-Ql 



2N r 



iVc 2 



C 2 (r)L(x,x) 



&f3a , 

(2.9) 
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(a) (b) 
Figure 3: Contractions between different Wilson lines; (a) ladder type and (b) crossing type. 



Figure 4: Typical diagram contributing to the Gaussian average of two Wilson lines, 
where we defined 

L(x,y) = g 4 J d 2 z±G (x ± -z ± )Go(y ± -z ± ) , (2.10) 
/v 2 — i r b ~ 

Q 2 s (b-,a-) = / dz-^(z-). (2.11) 



Here we remark that we will simply write Q 2 to denote Q 2 (+oo, — oo) in later discussions. 
It should be mentioned that L(x, x) does not depend on x in fact because of translational 
invariance, and thus we can write it as L(0, 0) equivalently. As we will argue later, however, 
L(x,y) generally suffers infrared singularity, and the expectation value ( |2.9| ) turns out to 
be negligible small. This observation intuitively corresponds to the fact that a single quark 
or gluon with non-trivial color charge would interact with color charge fluctuations inside 
the target at even far distance on the transverse plane. As a result of this long-ranged 
interaction (which should be cut off by either the target size or confining scale ~ ^qcd)' 
a quark or gluon is absorbed strongly in multiple scattering. 



2.2 Two-point function; (U(xx±)p iai U(x2±)p. 



2&2 I 



We shall next consider the two-point function of the Wilson lines. This is, in contrast to the 
one-point function, physically relevant if the projectile is a qq mesonic or gg glueball-like 
state. Not only the tadpole diagrams consisting of Fig. § (a) but also the ladder diagrams 
structured with the sub-diagram Fig. ||| (a) contribute to the Gaussian average (see Fig. § for 
typical example). Any diagram containing the crossing sub-diagram as shown in Fig. ^ (b) 
vanishes in the same way as the nesting and overlapping ones displayed in Fig. & We can 
compute the Gaussian average of two Wilson lines diagrammatically by combining all the 
tadpole and ladder subparts up. The Dyson equation we need to solve is as follows; 
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The horizontal lines represent (U(b ,a \x±)) involving all the tadpole contributions. We 
can literally convert the above graphical equation into the algebraic equation as 

(U (b~,a~ \xi±)/3 iai U (b~,a~ \x 2± )p 20i2 ) = 

= (U(b-,a-\x l± )) piai (U(b-,a-\x 2± ))p 2a2 - 



-g* / ^-(CTCft-.z-l^)) (^(t-.z-l^^t^^^MV) >< 

J a~ 

x J dz 2 ± G (x 1± - z ± )Go(x 2 ±- z±)(U(z~ , a~\x 1± ) ll0ll U(z~ , a~\x 2 ±) 120l2 ) 



U(b ,a \x 1± )U(b ,a \x 2± ) \ 6p iai 5p 2(X2 - L(x%, x 2 ) x 



x / dz-^(z-)t a pill t 



a 

(U(z- ,a \xi±)~ ilQl U(z ,a \x 2 ±).y 2a2 ) 



P2-12 ' 



■ (2-12) 



U(z ,a \x\±) U(z ,a \x 2 ±) 
We divide the both sides of this integral equation by U (b~ , a~ \x\±_) U (b~ , a~ \x 2 ±) to reach 
F(x\±, x 2 ±\b~ , a~)i3 lf 3 2 - ai a 2 = 

= S/hoH&faaz ~ L(xi,x 2 ) I dz n 2 {z )tp 1J1 t% 212 F(xi±,x 2 ±\z ,a ) 7l72; a ia2 , (2.13) 

J a~ 

where we defined 

(U(b~,a~\xi±)p iai U(b~,a~\x2±)p 2 a 2 ) 



F(x 1 _ L ,x 2± \b ,a ) / 3 l/ 3 2 -ai 



a 2 



U(b ,a \xi±)U(b ,a \x 2 ±) 



(2.14) 



The integral equation (2.13) has an identical structure as Eq. (|2,8|), so that we can solve it 
in the same way to find 



F(xi±,x 2 ±\b ,a )(3 1 (3 2 ;a 1 a 2 = exp 



-2Q 2 s (b~,a-' 



2N, 



' tf t 2 L(xi,x 2 ) 



m-i 



(2.15) 



/3i/32;aia2 



In the limit of a — ► —00 and b — > +00, as a result, we can express the Wilson line 
average in the following form; 



(U(x 1± ) /3iai U{x 2 ±) l 3 2 
= U{x 1± )U(x 2 ±)exp 



2N 

-IQl^^titlLixux-z) 



m - 1 



cxp 



exp 



Q 



2 2N C 
AT?-1 



2t\ t% L( Xl ,x 2 ) + tf L(x u a*) + tf L(x 2 ,x 2 



Pll32\OilCt2 



o2 



/3i/?2;cno2 



" Q - W=T + t2)2m 0) " t? * 5r(xi ' X2) 



(2.16) 



/3i/32;oi02 



where we made use of translational invariance to make a shift L(x, x) —* L(0, 0) and defined 
F(xi,x 2 ) = 2(L(0,0) — L(x\,x 2 )) which is free from infrared singularity. 
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Here, let us introduce "Hamiltonian" by 

(U(b~,a-\x 1± )p iai U(b~,a~\x 2 ±)(3 2a2 ) = exp[-(iJ + V)]^p 2 - aia2 (2.17) 
with the "free" part, 



and the "interaction" part, 



2/V 

H ° = ^iv^T (t? + tl)2m 0) ' (2 - 18) 



27V 

V = -Ql W ^l t i t 2 r ^2). (2.19) 



One can readily prove that [Hq, V] = meaning that THqT^ 1 and TVT~ l are to become 
diagonal simultaneously. Let us consider the decomposition of Hq into irreducible represen- 
tations in color space. We see that Hq is proportional to the second-order Casimir operator 
(tf + if) 2 . In general the irreducible representation is labeled by a set of non-negative in- 
tegers m with rank N c — 1, i.e., Dynkin coefficients. The associated second-order Casimir 
invariant is expressed as 



C 2 (m) = — 



iVc-l iVc-l 

m n )m n + 2 ^ l(N c - n)m n mi 

, n=l n>l 



(2.20) 



For example C 2 = 4/3 for the fundamental representation (triplet) of SU(iV c = 3) charac- 
terized by m = [1,0] and C 2 = 3 for the adjoint representation (octet) characterized by 
[1,1]. It is obvious from Eq. ( |2T20| ) that C 2 (to) is semi-positive and zero only when 



m 



to = 0, that is, a singlet. If the color structure of the Wilson line correlator is projected 
onto non-singlet states, Hq gives a large suppression factor, which can be seen from 



L(x,y) = g 4 J d 2 z ± Go(x ± -z ± )G (y ± -z ± ) 



(2vr)2 J (2vr)2 k 2 ± q\ 

4 /• dfcdfl _ a . JL | AicOB e l 
3 J (2vr) 2 /fc 3 

= Sf" 70 ^"^' (2 ' 21) 

where Jo(x) is the first-kind Bessel function. When k is small (strictly speaking, when 
k\x± — y_|_| is small), Jo(A;|a;_L — | ) ~ 1 — \(k\x± — y±\) 2 , hence, the momentum in- 
tegration in L(x,y) infraredly diverges. We introduce an infrared cutoff Aqqd to regu- 
larize the infrared singularity. Then, we have L(x,y) ~ V^qcd and F(x,y) ~ \x± — 
y 1 _\ 2 ln(\x± — y ± \AQCY)). Therefore, for small Aqcd, Hq which is proportional to L(0,0) 
should be much larger than V which is proportional to T(x,y). The non-singlet part in the 
color decomposition is thus accompanied by a large suppression factor, exp[— (2N C /(N 2 — 
1))C2(to)Q 2 /(47tAq CD )]. This is a physically reasonable result; the Wilson lines form a 
color singlet, and such neutral objects are free from long-ranged color interactions (except 
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for logarithmic singularity) which should be cut off by the confining scale Aqq D - Because 
Q s (which is ~ GeV order) is typically greater than Aqqd ~ fm _1 by one order of magni- 
tude at least, we do not have to concern non-singlet parts. From now on, accordingly, we 
consider only the singlet part of the color structure of the Wilson line product. 

2.3 n-point function; (U(x 1± )f 3iai U(x 2 ±)p 2a2 • • • U(x n± )f 3nCtn ) 

Now that we have understood the computational procedure, it is easy to generalize the 
formulae to the n-point case. The integral equation can be diagrammatically represented 
as 



The R.H.S. consists of all the possible permutations of one ladder bridging over two out 
of n Wilson lines. The horizontal lines are (U(b~ ,a~\x±)} with all the tadpole insertions. 
The corresponding equation is 

{U(b~ ,a~\x 1± )f3 iai U(b~ ,a~\x 2 ±)/3 2 a 2 ■ ■ ■ U(b~, a~\x nl _ )p n a n ) 

= Yl( u ( b '^ a ~\ XiA -))p iCH ~^ Yl / dz ~^ 2 { z ^) / dz 2 _ L G (x i± -z ± )G (x j± -z ± ) x 
i=l ' i>j Ja " J 

n 

x \\{U{b-,z-\ Xil _)) PkXk 5 Mll ■ ■ ■ Sx^^ t a Xn . 6 Xj+llj+1 ■ ■ ■ 

k=l 

x(U(z~ ,a~\x 1 ±) liai U{z~ ,a~\x 2 ±) l2a2 ■ ■ ■ U(a~, z~\x nl _)~ jn0ln ) . (2.22) 
The above integral equation takes a similar form to the case of two-point function in 



Eq. (2.12), so that we can find the solution in the same way as 

(U(x 1± ) (Siai U(x 2 ±)f3 2a2 ■ ■ ■ U(x n± ) PnCtn ) = exp[-(F + K)W- -fern •■•<*„ (2-23) 

with 

IN 

v = - Q2 ^w~i E ^ ^ r{x%u Xj±) ■ (2 - 25) 

c i>j 



These expressions are plain generalization of Eqs. (2.18) and ( |2.19| ). Again, Hq is propor- 



tional to the second-order Casimir operator, and after decomposing the color structure into 
irreducible representations, we can drop non-singlet parts. What we should do to simplify 
the result further is find the singlets out of the direct product of SU(iV c ) matrices which 
make Hq vanishing. For example, in case of the four-point function of Wilson lines in the 
adjoint representation with N c = 3 (i.e. four gluon propagation through a dense target), 
there are eight independent singlets out of 8 ® 8 (g) 8 (g> 8. Thus, the Wilson line correlator 
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can be expressed to be an 8 x 8 matrix of e - ^ in the basis of eight singlets. We will face 
with concrete calculations later. In most cases of our interest in physics problems, all we 
need to know is expressed in terms of its eigenvalues and eigenstates. 

3. Examples 

We will elaborate several concrete examples relevant to physical processes. We will identify 
the singlet basis to compute the eigenvalue of the color matrix. 

3.1 (U{x 1 ±) / 3 iai U*(x 2 ±)p 2 a2 ) in the fundamental representation 

In our formulae fl2.18| ) and fl2.19| ) we set tf = Tp and t% = — Tp* where Tp's are the 



SU(JV C ) generator in the fundamental representation. This expectation value appears in 
one color dipole or qq scattering off a dense target. In this case the number of singlet 
is only one; N c <S) N* = 1 © N% — 1. If we denote the singlet state as \s) then we have 
CK2 1 — S aiCt2 / \J Nc with a proper normalization. The projected element of the necessary 
part is 

(s\v\s) = Ql^i r F a , iai T ¥ f 2a2 r( Xl , X2 ) 5 ^Mi = qI r( XuX2 ) . (3.1) 

Consequently, the singlet part of the two-point function of Wilson lines in the fundamental 
representation is 

(U(x 1± )( Siai U*(x 2 ±)p 2 a 2 ) = exp[-Qlr(x 1 ,x 2 j\ 5p ll32 5 aia2 /N c , (3.2) 

which means that only the closed color dipole survives under the average over random color 
distribution inside a dense gluon medium. 

3.2 {U(xi±)p iai U(x 2 ±)/3 20t2 ) m the general representation 

It is possible to extend the previous argument formally. We can reach a more general ex- 
pression in arbitrary representation. The direct product of two irreducible representations, 
r\ and r 2 , makes one singlet when r 2 is the star representation to n, that is, r± = r\ = r, 
leading to the generators, t% = —tf*. In this general case we have (aia 2 \s) = S aia2 /y/d^, 
where d r is the dimension of r representation, which is given by the Dynkin coefficients m 
charactering the irreducible representation r as 

Af c / i-i \ 
d r (m) = H U + E^7 • M 

i<j \ n=i J / 
We can express tft^ using the Casimir invariant; 

«?<S = \ [(<! + t2? ~ ~ W) 2 ] = -\[C2(n) + C 2 (r 2 )] = -C 2 (r) , (3.4) 

IN 

(s\V\s) = Qlj^ C 2 {r) r( Xl ,x 2 ) . (3.5) 



leading to 



After all, the generalization of Eq. (|3.2| ) takes an expression of 

T 2N 

<^(a;i±)/3 iai ^(«2±)/3 2 Q 2 ) = exp -Ql n2 ^ - C 2 (r)r(xx,x 2 ) 



Splat 5p 2 a z /d r . (3.6) 
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3.3 (U(xi±) ai p 1 U(x2±)a 2 i32 ' " U(xN c ±) aN p N ) in the fundamental representation 

This expectation value is relevant to the scattering amplitude between a baryon consisting 
of N c valence quarks and dense gluon matter inside a heavy hadron. One can find an 
example of the baryon expectation value in the context of the Odderon physics [fl4|, 15], 
though the Gaussian average can only describe the Pomeron part. 

The singlet state with one baryon is (a\ ■ ■ ■ a n \s) = e aia2 ... aN /y/N c \ where e aia2 ... aN is 
the antisymmetric tensor with the definition ei2—jv c = 1- The matrix element (s|Tf"Tf"|s) 
does not depend on the indices, i and j, that is, 

(s\T F ? T^\s) = "^f^ (N c ~ 2)1 = -^±1 , (3.7) 

that yields 

(s\V\s) = Ql — Yl r(x i± ,x j± ) . (3.8) 
The baryon expectation value is thus, 

(U(x 1± )p iai U{x 2 ±)p 2 a 2 ■ ■ ■ U{x Nc± )p NcaNc ) = 

-Q* J[~L E ea ^ N J^ . (3.9) 

When N c gets large, the exponential factor is decreasing in contrast to the meson scattering 
in Eq. ( |3.2| ) that stays unsuppressed for large N c . This is a manifestation of the fact that 
baryons would not live as they are in the large- iV c limit but mesons would. 

3.4 (U(xi_ L )p lQl U(x 2 ±)p 2a2 U(x 3 _ L )p :jas ) in the adjoint representation 

So far, the number of singlet is only one, and the next step we are heading for is to 
treat the case with multiple singlets. The simplest and still non-trivial is the three gluon 
propagation, in which — 1) (8) (N% — 1) ® (N% — 1) includes two singlets composed of 
f abc = -2itr([r|,T|]TX) and d abc = 2tr({T^, T%}T£), where T£'s represent the SU(7V C ) 
generator in the adjoint representation, as 

1 



exp 



(aia 2 a 3 \si) = f ai 



OL2Q-3, 



- i)jvc 



(aia 2 a 3 |s 2 ) = d aia2a3 J ^ _ _ fj ( 3 - 10 ) 

with a proper normalization. Using the Jacobi identities, fabefcde + fadefbce + facefdbe = 
and fabedcde + fadedbce + faceddbe = we can calculate each matrix element in singlet space 
as 

1 N 

( Sl \T A tT A ]\ Sl ) = -f a 

(si\T A f Txj\a 2 ) = -fcx 1 ap 1 fc l2 af3 2 fa 1 a 2 l3 3 dp l p 2 p 3 — — , - = 0, 

(s 2 \T A f T A °j\s 2 ) = -fa 1 a(3 1 fa 2 ap 2 d aia2 p 3 dp 1 p 2 p 3 _ 4) (^2 _ 1) = ~ ~Y ' ( 3-11 ) 
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which are independent of the indices i and j, and thus the matrix structure is simply 
proportional to unity. The projected V is then, 



N 2 



x 



i>3 



1 0' 

1, 



(3.12) 



At last, we get the result, 

(U {Xi±)i3 iai U(x 2 ±)f3 2 a 2 U {x 3 ±)fr 



3 a 3 / 



N C (N? - 1) 



exp 



■Q 



2 ^ 



( N ' 



AT2 



-d, 



. (3.13) 



Roughly speaking, three gluons can make two kinds of color singlet glueballs and they 
do not mix together due to different symmetry. The first and second terms in the curly 
parenthesis are contributions from those glueball states respectively. 

3.5 (U(xi±) j 3 iai U*(x2±)f3 2 a 2 U( K x 3 ±)f3 sas U*(x 4: ±)f3 4a4 } in the fundamental representa- 
tion 

The next non-trivial example is the four-point Wilson lines in the fundamental representa- 
tion. This problem reduces to the irreducible decomposition of N c ® N* ® N c ® N*, which 
contains two normalized singlets; 



(aia2<y-3<y-4\si) 
{ot\oi2aza4\s2) 



1 



^a±a 2 "i 



1 



leading to the projected matrix elements of 
2N C 



a 2 as 



0304 I j 



(3.14) 



V 



Q* j^rzi ( - t f? T F2* r(x 1± ,x 2 ±) + t f i T F a 3 r( Xl± ,x 3± ) - t f ? T F a * r(x 1± ,x 41 _) 



-Tpf T F a 3 r(x 2± ,x 3± ) + T F 2* T F ? r(x 2± ,x 4± ) - T F a 3 T F f r(x 3± , X41) 
given by 

/ 



(3.15) 



V n V 12 
V21 V22 , 



n 2 2N C 
N} - 1 



o 



(3.16) 



where we defined with slight modification from Ref. 

2a = r(x 1± , x 41 _) + r(x 2 ±, x 3± ) 

2(3 = r(xi±,x 3 ±) + r(x 2 ±,x 4± ) 
27 = r(x 1± ,x 2 ±) + r(x 3± ,x A1 _) 



(3.17) 
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which are free from infrared singularity. In order to calculate the matrix element of 
exp(— V), we need to diagonalize V whose eigenvalues are given as the solution of the 
characteristic equation det[A — V] = 0, which gives us two eigenvalues, 



1 



\ ± = -(tTV±<p) 



(3.18) 



with 



. IN 2 4 

<p = v / (trl/) 2 -4dety = Q 2 sW ±j J (a - 7 ) 2 + ~ a)(/3 ~ 7) (3-19) 



in our notation. The eigenstates are 



u+ 



(cos 6 s 
I sin0 




f-sin0\ 

l COS0 J 









(3.20) 



with 



tan ( 



Vn -V 2 2 + <p 

Using the unitary matrix T = (u + u_) we can diagonalize exp(— V) to have 



(3.21) 



e" 







\ e5*y 
tr y | cosh I9 9 — (cos 2 — sin 2 9) sinh iy? 



_ i 
e 2 



e -ftrV 



cosh |</? 



-2 cos sin sinh ±<p 
sinh 1^3 / Vu - V22 2V12 



— 2 cos 9 sin sinh \ ip 
cosh + (cos 2 9 — sin 2 0) sinh |<£ 



2Fi 2 -Vu + V22 



(3.22) 



We can find one direct application of this example in the quark production from the 
CGC background in the p-A collision ||, in which the necessary quantity is 

(tr{Z7(aJi_L) T§ U^x 2 ±) U(x 31 _) T§ U\x 4± )}) = 

= r Fj3 1 /3 2 r Fj33/34^a2a3^4ai(^( S l±)ai/3i^*(^2±)a 2 /3 2 ^(^3±)Q 3 /33 C/ *( a; 4±)a4/34) ■ (3-23) 

The color structure in ^FAftj^Fflj&^oa^ai is to be expressed in terms of the singlet 
basis of Eq. ( p. 14 ) as follows; 



rp a rp a r c 
l¥ 0i02 l¥ 0Z04 "2"3 °a 4 ai 



{Pi(32P3^\s2)^V N c ~ l(s2\aia 2 a 3 a 4: ) + {si\aia 2 a 3 a^)\ . (3.24) 



- 12 - 



Thus we immediately conclude 



(ti{U(x 1± ) Tp U\x 2± ) U{x 31 _) T§ C^(x4l)}) 



(s 2 |e v |s 2 ) + 



1 



( s l| e \ s 2, 



N?-l 



•e-5 try 



cosh \ip 



sinh iy; 



22 



12 



1 r 



Q 2 2N 2 
cosh±^ + -± jpZi ^-u) sinhi^j 



. (3.25) 



In order to compare this result to Eq. (86) in Ref. f|, we note that Q 2 = (N 2 - l)/(AN c )/j,\ 
and ip = \N c ^x 2 A \a — 7|vA if written with the notation of Ref. ||. Then, our result dj|D 
turns out to agree exactly with Ref. ||. 

3.6 (U{x 1± ) /3iai U(x2±)p 2a2 U(x 31 _)fy aa JJ(x 4:± ) l3ia4 ) in the adjoint representation 
This quantity is, literally speaking, the four-gluon scattering amplitude. Also, in the 



calculation of two-gluon production from the CGC background [11], we have to evaluate 
tr{U T^U^U T£U^} in the adjoint representation. There are eight singlets [16] given by 
combinations of the following bases; 



^ai«27 ^03«47 ^aia37^«2a47 
^01027 ^03047 ^ctict37 J02Q47 ^aict47 ^02037 



(3.26) 



These states are not orthogonal, and it is very hard to obtain the eigenvalues and eigenstates 
of an 8 x 8 matrix V to compute e~ v . In the next section, we will develop an approximation 
to simplify the calculation. 



4. Large-A^ limit 

It is possible to construct any higher-dimensional representation from the direct product 
of the fundamental (and anti-fundamental) representation. A well-known example is the 
adjoint representation whose matrix representation can be given in terms of the fundamen- 
tal Wilson lines as 2tr[U(x±)t^U^(x±)t a ] where a and [5 run from one to N 2 — 1. Here 
we shall consider the arbitrary product of Vs in a general way in the large-iV c limit. As 
usual, then, among color singlets, the baryon operator is dropped and only the meson-type 
operators remain non- vanishing, as we have seen before. That is, we shall focus on the 
singlets out of N c <g> N* <g> ■ ■ • <g> N c <g> N*. 
The Gaussian average of our interest is 

(flUix^U*^)-^) = exp[-(F + ^W...^^,..^ , (4.1) 
i=l 
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where we can write 



//<* - Q s 2 t^tL(0,0) 



V 



-Ql 



N*-l 
2N< 



i=i 



(4.2) 



-£r F ?T F f r(x i± ,y j± ) , 
t,i=i J 

before taking the limit of large iV c . The second-order Casimir operator in H$ is 

2 



(4.3) 



£(T F ?-T F f 



i=l 



J2(T F f+T F f 2 )+2 y }2(T F *T F «+T F fT F f)-2 ^ T F ?T F f . (4.4) 



i=l 



The first term of R.H.S. in Eq. ( |4,4j ) is just the Casimir operator and thus proportional to 
a unit matrix, 



Y,(T F f + T F f 2 ) = 2n- 



i=i 



2N r 



in the large- N c limit. Next, let us check that the second term of R.H.S. in Eq. 
to only 0(1) contributions. A well-known formula reads 



rp a rp a 



Pj a 3 



1 



N, 



if/% 



(4.5) 
leads 

(4.6) 



This cannot make a loop with any singlet state; a loop formed with ati-Pj and Oij-/?i 
connected is not a singlet with respect to q,'s or (3^s which should be accompanied by a 
huge suppression factor in the infrared sector. This is why the second term is only negligible 
in the large- N c limit. The last term of R.H.S. in Eq. (4.4) has 0{N C ) contributions as seen 
from 



that can make a loop like 



rp CI rp a* 



2 VA 



rp CI rp a* ? 



2VA- 



2 VA 



(4.8) 



From Eqs. (4.5) and (|4.8|) , we can find the color bases on which Hq vanishes. There are n! 
relevant singlets given by all the permutations of 5 aiSll • • • S a „a n , i.e. 



(ai • • • a n ; «i • • • a n |s p ) 



(4.9) 



where (p%, . . . ,p n ) is the permutation of (1, ... , n) labeled by p which runs from one to n\. 
Here we remark that we can relax large N c for Eq. ( [4,9|) being a singlet, though we derived 
it in the large- N c limit. 
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It is easy to evaluate V in this basis because the color matrix structure of V is quite 
similar to Hq at large N c . So that, we have 



V\s p 



2Q| 

N 



J2TtfT F fr(x i± ,y j± )\s p ) -» Q 2 s ^2r(x t± ,y pi± )\s p ). (4.10) 



1=1 



where we used Eq. (|P]). We finally arrive at the expression of the Gaussian averaged 
Wilson loops in the large- N c limit as follows; 



flU(x i± )p iai U*(y 



i=l 



1 

iv: 



(4.11) 

We shall apply our formula (4.11) for evaluation of the expectation value of the dipole 
operator defined by 



D(x±,y ± ) = —tr[U(x ± )U^(y ± )] = —6 p p 6 a5l Uix^pJJ*^^^ . 



Nr. 



The scattering amplitude of the light projectile with n color dipoles is written as 



Y[D(x i± ,y i± )\ = exp 



8=1 



t=i 



(4.12) 



(4.13) 



It should be mentioned that, when y i _ L — > Xi±, the exponential factor becomes one with 
r(xi±,y i _ L — ► a^jj — ► 0, meaning color transparency. 

We shall next compute the n-point Wilson line correlator in the adjoint representation, 
which is as easy as 



i=l 



Y[U(x nX )b iai 

n 



i=l 



n 



(4.14) 



Finally, we make a comment on the calculation of two-gluon production from the CGC 
background. In such a case, unfortunately, the leading-A^ order of tt{UT?U<UTfU'} is 
just vanishing and the sub-leading order is necessary. We will report details elsewhere p| . 

5. Summary 

We have derived the general formula to compute the correlation function of Wilson lines 
in the random distribution of color source, i.e. in the McLerran-Venugopalan model. We 
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emphasize that our technique is to be applicable whenever the CGC weight function is 
approximated as a Gaussian. Our formula would be quite useful in calculations not only 
of the scattering amplitude but also of the particle production from the CGC background. 
The correlation function or the scattering amplitude is strongly suppressed if the color 
non-singlet irreducible representation of the n-particle initial or final state is involved. 
Hence, we only have to consider the singlet part of the color structure associated with the 
Wilson line product in order to evaluate the correlation function. After all, the problem of 
evaluation of the correlation function is simply reduced to diagonalization of a color matrix 
using singlet bases. We have explicitly written the two-point function down in the general 
representation, the three-point function in the fundamental representation corresponding 
to baryon's scattering off the CGC background, in the adjoint representation as well, and 
the four-point function in the fundamental and anti-fundamental representation which is 
related to the qq production from the CGC background. 

The larger number of Wilson lines, n, is involved in the Gaussian average, the more 
difficult it is to find the singlets and to diagonalize the color matrix. As a matter of fact, the 
number of the singlet states increases exponentially with increasing n, although our method 
is powerful enough to implement also in numerical computations in such an intricate case. 
Instead of that, we took advantage of simplifying the expressions in the large- N c limit. We 
have derived the explicit formula with arbitrary number of Wilson lines in the fundamental 
representation at large N c , from which, in principle, any representation can be constructed. 
For the phenomenological application, we plan to address our calculation of the two-gluon 
production from the CGC background in another publication jHJ where we will discuss 
the forward-backward rapidity correlation with respect to the hadron multiplicity in the 
collision. 
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